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Abstract 

We observe that the saturation model of deep inelastic scattering, which suc- 
cessfully describes inclusive and diffractive data at small x, predicts a geometric 
scaling of the total 7*p cross section in the region of small Bjorken variable x. The 
geometric scaling in this case means that the cross section is a function of only one 
dimensionless variable r = Q'^Rq{x), where the function Ro{x) (called saturation 
radius) decreases with decreasing x. We show that the experimental data from 
HERA in the region x < 0.01 confirm the expectations of this scaling over a very 
broad region of Q^. We suggest that the geometric scaling is more general than the 
saturation model. 
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It has recently been observed that the ep deep inelastic scattering (DIS) data at low 
X m], 0] can be very economically described with the help of the saturation model 0]. 
In this model the QCD dipole picture of the interaction between the transversely (T) or 
longitudinally (L) polarized virtual photon 7* (emitted by the incident electron) and the 
proton p was adopted. In this case the total 7*p cross sections are given by |§, ||, H 

aT,L{^,Q^) = ldh£dz\<ilT,L{'r,z,Q')\'Hr,x), (1) 

where \I/t,l is the wave function for the splitting of the virtual photon into a qq pair 
(dipole), and a is the imaginary part of the forward scattering amplitude of the qq dipole 
on the proton, called the dipole cross section, which describes the interaction of the 
dipole with the proton. In addition, r is the transverse separation of the quarks in the qq 
pair, and z is the light-cone momentum fraction of the photon carried by the quark (or 
antiquark). As usual, —Q^ is the photon virtuality and x is the Bjorken variable. The 
cross section (|I|) has a clear physical interpretation in the proton rest frame in which the 
qq pair formation and a subsequent interaction with the proton are clearly separated in 
time. 

Let us recall that the standard DIS proton structure functions are related to a^^L by 

FtA^, Q^) = GT^X, Q2) , (2) 

47r ttem 

and F2 = Ft + Fl. The wave function of the virtual photon is given by the following 
equations: 

+ m} K'.iQfr)} , (3) 

where the sum is performed over quarks with flavour /, charge Cf and mass ruf, and 

= z{l - z)Q^ + m) . (4) 

The functions Kq i are the Bessel-Mc Donald functions. 

The main assumption of the saturation model concerns the saturation property of the 
dipole cross section which is incorporated in the approach of ref. 0| as below: 

a(x,r) = ao,(^). (5) 

The function Rq{x) with the dimension of length, called saturation radius, decreases with 
decreasing x, while the normalization ctq is independent of x. When f = r / Rq{x) —>■ 00 the 
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function g{f) saturates to 1, so that r) — > (Tq. In the reahzation of the saturation 
model 

g{f) = 1 — exp (— . (6) 

The fact that the dipole cross section (j^) is limited by the energy independent cross section 
(Jo may be regarded as a unitarity bound. This reflects the fact that the small x increase 
of DIS structure functions generated by pure DGLAP or linear BFKL evolution has to 
be tamed by unitarization effects |7|]- [p!^ . 

The characteristic feature of eq. (|^) is its "geometric scaling", i.e. cr{x, r) only depends 
on the dimensionless ratio r/Ro{x), and its energy dependence is entirely driven by the 
saturation radius Ro{x). The scaling property of (H) with ctq independent of x resembles 
geometric scaling of hadron-hadron scattering |^D|] . In the latter case the relevant quantity 
is the scattering amplitude in the impact parameter representation G{b'^,s), where b is 
the impact parameter and s is the CM energy squared. The geometric scaling in this case 
corresponds to the assumption that 

G{b\s)-^Gi(3), (7) 

where f3 = b'^/R^{s) with R{s) corresponding to the interaction radius which increases 
with increasing energy. The analogy between scaling exhibited in hadron-hadron collisions 
and in deep inelastic scattering should not be taken too literally. For example, the two 
radii have different physical interpretation, and moreover, they show completely different 
energy dependence since the saturation radius Ro{x) decreases with increasing energy (for 
X ~ Q'^/s 0). 

The assumption about the scaling property of the dipole cross section (|^) has profound 
consequences for the measured 7*p cross section a^*p = ar + c^l- If we neglect the quark 
masses in the photon wave functions (|^) we can rescale the dipole size r — > r / Rq{x) in 
eq. (|1|) such that the integration variables are dimensionless. Thus, after the integration 
a^*p becomes a function of only one dimensionless variable r = Q'^R'^{x), instead of x and 

separately, 

where the scale for a^^p is provided by the dipole cross section normalization o"o- The 
non-zero light quark mass does not lead to a significant breaking of the scaling (|^) . Let us 
emphasize again that the new scaling for cr^*p, valid in the low x region, is obtained due 
to the assumption that (H) depends on r and x through the dimensionless combination 
r/Ro{x). Following the discussion in it is easy to show that with the realization (||) 
we smoothly change the behaviour of (||), 

0"7*p ~ O"0 ^ 0"7*p ~ O'd/t , (9) 

(modulo logarithmic modifications in r) when r changes from small to large values, respec- 
tively. The aim of this paper is to demonstrate that the DIS data do indeed approximately 
exhibit the geometric scaling (W) with the property (0). 
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Let us discuss at first the parameterization of tlie saturation radius Ro{x). In ref. ||^ 
tlie saturation radius form was postulated as follows 

Mx) = TT -) ' (10) 



Qo Vxo/ 

where Qo = 1 GeV. The parameters Xq, A > together with the dipole cross section 
normalization o"o were fitted to all inclusive DIS data with x < 0.01. A very good de- 
scription of data was obtained down to the region Q"^ < 1 GeV^. The saturation model 
was in particular able to describe the transition from the region of DIS to the region of 
low values of Q^. Also DIS diffractive data were described without additional tuning of 
the parameters. For a recent related analyses see and also . 



Below we show that Ro{x) can also be determined differently in a less model-dependent 
way. To this aim let us observe that after suitable extension of the saturation model to 
the low region including the photoproduction limit Q^ = 0, the x dependence of 
the saturation radius Ro{x) can be correlated with the energy dependence of the total 
photoproduction cross section a^p. In order to extend the saturation model to the region 
of low Q^ we replace, following ref. [Q, x by a: defined as: 

^ = ^(l + ^) (11) 

in the argument of Ro{x) in eq. (^) and keep ruf ^ 0. The variable x is related to the 
total energy W by 

We note that the saturation model based on eqs. (|l])-(|D can now be extended down to 
the region Q^ = 0. The photoproduction cross section is given by equations with 
Q^ = 0, = Tn'j and with x replaced by a; = Arnj/W"^. The dominant contribution to 
the photoproduction cross section comes from the integration region l/m'j ^ ^ -Ro(^) 
in the corresponding integral on the right hand side in eq. (|lD. In this region we can 
set m'jKf[mfr) ~ and a{x,r) ~ ctq. This gives the following relation between 

photoproduction cross section and the saturation radius 

The parameter ctq is related to the overall normalization of the dipole cross section ctq by 
(Jo = (2aem,/37r)cTo. From equation we finally obtain the following prescription for 
the saturation radius 

Rlix) = ^ exp (-^) . (14) 

For a.yp we take the Donnachie-Landshoff parameterization ^3 

a^p = a , (15) 
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where we set = 140 MeV (following 0) in eqs. (|T^ and (|T^). Using results of the fit 
presented in |^ we find 

2 \ 0.08 

(16) 



/ 4m2 V 



For (To we set 23 fib to obtain a good description of data. 

Let us now confront the implications of geometric scaling (^ with experimental data 
on deep inelastic scattering at low x. In Fig. 1 we show experimental data |l[] on the total 
cross section cr^.p plotted versus scaling variable r = Q'^Rl{x), with Ro{x) obtained from 
Eq. (|I4[). We include all available data for x < 0.01 in the range of values between 
0.045 GeV2 and 450 GeV^. 

We see that to the data exhibit geometric scaling over a very broad region of Q^. 
We can also clearly see the change of shape of the dependence of a^p on r from the 
approximate 1/r dependence at large r to the less steep dependence at small r. Thus the 
asymptotic relations (|^) are also to a good approximation confirmed. The approximate 
asymptotic 1/r dependence reflects the fact that the cross section a^*p scales as 1/Q^ 
(modulo logarithmic corrections) and its energy dependence is governed by 1/Rl{x). Less 
steep dependence corresponds to the fact that at small values of r the total cross section 
grows weaker with energy than 1/Rq{x) due to saturation of the dipole cross section, see 
eq. d). 

We also found a symmetry between the regions of large and small r for the function 
A/rcr^*p, which is illustrated in Fig. 2. For the asymptotic values of r this is a manifestation 
of the relations (^). It is remarkable that Fig. 2 seems to indicate the presence of symmetry 
of y/Ta^*p with respect to the transformation 

r < — > 1/r (17) 

in the whole region of r. 

We have also tried the power law parameterization for the radius, 

Rl{x)^x\ (18) 



where 0.3 < A < 0.4, in particular the original form proposed in 0] (see eq. (p!0[) ), and 
found that the data also exhibit the geometric scaling with this choice of parameterization. 
The approximate 1/r dependence at large r corresponds to the behaviour of the 
proton structure function F2 at large Q^. 

In the photoproduction case, the parameterization (pIS]) combined with relation ([T3|) 
would correspond to the logarithmic dependence on energy of the photoproduction cross 
section, i.e. 

a^p ~ ln(a;) (19) 



We do therefore find that both prescriptions for a^p eq. (|T5]) and eq. (|19D give numerically 
similar results. 
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In Fig. 3 we show contours corresponding to different values of variable r in {x, Q"^) 
plane together with experimental points for these values of r. Geometric scaling means 
that o"^*p is constant along each line. To be precise for each value of r we plot experimental 
points within the bin (ln(r) ±5) with 5 = 0.1. We see from this figure that for each r there 
are several experimental points for which x varies as much as two orders of magnitude 
and changes by a factor of 4. Despite of that, all points along each line in Fig. 3 are 
transformed to a narrow spread of points for a particular value of r in Figs. 1 and 2, thus 
exhibiting geometric scaling. 

In order to show that the geometric scaling is confined only to the small x region we 
show a similar plot as in Fig. 1, but for the experimental data with x > 0.01 [Q, |^. It 
is evident that the geometric scaling is significantly violated in the region of large x, see 
Fig. 4. 

We would like to emphasize that the (approximate) geometric scaling should predom- 
inantly be regarded as a remarkable regularity of DIS experimental data at low x. 

Although this scaling has been inspired by the saturation model of ref. [Q, its sig- 
nificance goes presumably far beyond this model. In its essence the new scaling is a 
manifestation of the presence of an internal scale (saturation scale) characterizing dense 
partonic systems, Qs{x) ~ 1/Rq{x). This scale emerges from a pioneering work of 0, 
which was subsequently analyzed and generalized in |]Ty|. In the analysis |T^, and 



more recently in |jT8| , the scaling properties similar to those postulated in (^ were found. 



An independent formulation of the small x processes, gives the same overall picture 
with the saturation scale. At a deeper level, the geometric scaling for small-x processes 
may reflect self similarity or conformal symmetry of the underlying dynamics. More 
detailed studies are under way, see |T6|-p9[]. 

To sum up we have shown that the experimental data on deep inelastic ep scattering 
at low X exhibit geometric scaling, i.e. the total cross section a^*p{x, Q^) is the function 
of only one dimensionless variable r = Q'^Rl{x). This regularity was found to hold over 
the very broad range of from 0.045 GeV^ to 450 GeV^. It would be interesting to 
understand in detail a possible dynamical origin of this simple regularity. 



Acknowledgments 

We thank Jochen Bartels and Genya Levin for useful comments and Larry McLerran for 
useful correspondance. This research has been partially supported by the EU Framework 
TMR programme, contract FMRX-CT98-0194, Deutsche Forschungsgemeinschaft and the 
Pohsh Commitee for Scientific Research grants Nos. KBN 2P03B 120 19, 2P03B 051 19. 



6 



References 



[1] HI Collaboration (S. Aid et al), Nucl. Phys. B 470 (1996) 3; ZEUS Collaboration 
(M. Derrick et al.), Z. Phys. C 72 (1996) 399; HI Collaboration (C. Adloff et 
al.), Nucl. Phys. B 497 (1997) 3; ZEUS Collaboration (J. Breitweg et al.), Phys. 
Lett. B 407 (1997) 432; ZEUS Collaboration (J. Breitweg et al.), DESY-00-071, 
|hep-ex/00050l4 

[2] BCDMS collaboration (A.C. Benvenuti et al), Phys. Lett. B 223 (1989) 485; NMC 
collaboration ( M. Arneodo et al.), Phys. Lett. B 364 (1995) 107; E665 collaboration 
( M.R. Adams et al.), Phys. Rev. D 54 (1996) 3006; SLAC collaboration, Phys. Lett. 
B 282 (1992) 475; EMC collaboration (Aubert et al.), Nucl. Phys. B 259 (1985) 
189. 

[3] K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D59 (1999) 014017; Phys. Rev. D60 
(1999) 114023. 

[4] J.D. Bjorken, J.B. Kogut and D.E. Sopper, Phys. Rev. D3 (1971) 1382. 

[5] N.N. Nikolaev and B.C. Zakharov, Z. Phys. C49 (1991) 607; Z. Phys C53 (1992) 
331; Z. Phys. C64 (1994) 651; JETP 78 (1994) 598. 

[6] A. H. Mueller, Nucl. Phys. B415 (1994) 373; A. H. Mueller and B. Patel, Nucl. Phys. 
B425 (1994) 471; A. H. Mueller, Nucl. Phys. B437 (1995) 107. 

[7] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1. 

[8] A.H. Mueller and Jian-wei Qiu, Nucl. Phys. B268 (1986) 427 ; A.H. Mueller, Nucl. 
Phys. B335 (1990) 115; Yu.A. Kovchegov, A.H. Mueller and S. Wallon, Nucl. Phys. 
B507 (1997) 367. A.H. Mueller, Eur. Phys. J. Al (1998) 19; Nucl.Phys. A654 (1999) 
37c; Nucl.Phys. B558 (1999) 285. 

[9] J. C. Collins and J. Kwieciriski, Nucl. Phys. B335 (1990) 89. 

[10] J. Bartels, G.A. Schuler and J. Bliimlein, Z. Phys. C50 (1991) 91; Nucl. Phys. Proc. 
Suppl. 18 C (1991) 147; J. Bartels, Phys. Lett. B298 (1993) 204; Z. Phys. C60 

(1993) 471; Z. Phys. C62 (1994) 425; J. Bartels and M. Wiisthoff, Z. Phys. C66 
(1995) 157; J. Bartels and C. Ewerz, JHEP 9909 (1999) 026. 

[11] J. Bartels and G. Levin, Nucl. Phys. B387 (1992) 617. 

[12] L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 2233, Phys. Rev. D49 

(1994) 3352, Phys. Rev. D50 (1994) 2225; A. Kovner, L. McLerran and H. Weigert, 
Phys. Rev. D52 (1995) 6231, Phys. Rev. D52 (1995) 3809; R. Venugopalan, Acta. 
Phys. Polon. B30 (1999) 3731. 



7 



[13] A. Biaks and R. Peschanski, Phys. Lett. B355 (1995) 301; Phys. Lett. B378 (1996) 
302; Phys. Lett. B387 (1996) 405; A. Bialas Acta Phys. Polon. B28 (1997) 1239; 
A. Bialas and W. Czyz, Acta Phys. Polon. B29 (1998) 2095; A. Bialas, H. Navelet 
and R. Peschanski, Phys. Rev. D57 (1998) 6585; Phys. Lett. B427 (1998) 147; 
|hep-ph/000924£ . 



[14] G.P. Salam, Nucl. Phys. B449 (1995) 589; Nucl. Phys. B461 (1996) 512; Comput. 
Phys. Commun. 105 (1997) 62; A.H. Mueller and G.P. Salam, Nucl. Phys. B475 
(1996) 293. 

[15] E. Gotsman, E. M. Levin and U. Maor, Nucl. Phys. B464 (1996) 251; Nucl. Phys. 
B493 (1997) 354; Phys. Lett. B245 (1998) 369; Eur. Phys. J. C5 (1998) 303; E. 
Gotsman, E. M. Levin, U. Maor and E. Naftali, Nucl. Phys. B539 (1999) 535; A. L. 
Ayala Filho, M. B. Gay Ducati and E. M. Levin, Eur. Phys. J. C8 (1999) 115; E. 
Levin and U. Maor, |hep-ph/0009217[ 

[16] la. Balitsky, Nucl.Phys. B463 (1996) 99. 

[17] J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D55 (1997) 
5414; J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D59 (1999) 014014; 
Phys. Rev. D59 (1999) 014015; Phys. Rev. D59 (1999) 034007; Erratum-ibid. 
D59 (1999) 099903; A. Kovner, J.Guilherme Milhano and H. Weigert, OUTP-00- 
10P,NORDITA-2000-14-HE, |hep-ph/0004014l; H. Weigert, NORDITA-2000-34-HE, 



liep-ph/0004044| . 



[18] Y.V. Kovchegov, Phys. Rev. D60 (1999) 034008; Phys. Rev. D61 (2000) 074018; G. 
Levin and K. Tuchin, Nucl. Phys. B537 (2000) 833; M.A. Braun, Eur. Phys. J. C16 
(2000) 337. 

[19] Y.V. Kovchegov and L. McLerran, Phys. Rev. D60 (1999) 054025; Erratum-ibid. 
D62 (2000) 019901; Y.V. Kovchegov and G. Levin, Nucl. Phys. B577 (2000) 221; 

[20] J. Bias de Deus, Nucl. Phys. B 59 (1973) 231; A. J. Buras, J. Bias de Deus, Nucl. 
Phys. B 71 (1974) 481; J. Bias de Beus, P. KroU, J. Phys. G 9 (1983) L81; J. Bias 
de Beus, Acta Phys. Polon. B 6 (1975) 613. 

[21] W. Buchmiiller, T. Gehrmann and A. Hebecker Nucl. Phys. B537 (1999) 477; G. 
R. Forshaw, G. Kerley and G. Shaw, Phys. Rev D60 (1999) 074012, Nucl. Phys. 
A675 (2000) 80; E. Gotsman, E. Levin, U. Maor and E. Naftah, Eur. Phys. J. CIO 
(1999) 689; M. McBermott, L. Frankfurt, V. Guzey and M. Strikman, Eur. Phys. 
J. C16 (2000) 641; T. G. Cvetic, B. Schildknecht and A. Shoshi, Acta Phys. Polon. 
B 30 (1999) 3265; A. Capella, E. G. Ferreiro, A. B. Kaidalov and C. A. Salgado, 
|hep-ph/0006235 . 



[22] A. Bonnachie and P.V. Landshoff, Phys. Lett. B437 (1998) 408. 
[23] A. Bonnachie and P.V. Landshoff, Phys. Lett. B296 (1992) 227. 



8 




Figure 1: Experimental data on a^»p from the region x < 0.01 plotted versus the scaling 
variable r = Q'^Rl{x). 
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Figure 2: \/rcr^*p plotted versus the scaling variable r 
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Figure 3: The lines corresponding to different values of scaling variable r (continues 
curves) in the (x, Q^)-plane. The points correspond to available experimental data located 
within the bins ln(r) ± 6 {S — 0.1) for each value of r. The numbers correspond to the 
value of r for each curve. 
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Figure 4: Experimental data on a^*p from the region x > 0.01 plotted versus the scaling 
variable r = Q'^Rl{x). 
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